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We evaluate a number of new instanton-induced phenomena in QCD, starting with static dipole- 
dipole potentials, and proceeding to quark-quark and dipole-dipole scattering at high energy. We 
use a non-perturbative formulation of the scattering amplitude in terms of a correlator of two 
Wilson-lines (quarks) or Wilson-loops (dipoles) and analyze the Euclidean amplitudes with both 
perturbative gluons and instantons. The results are analytically continued to Minkowski geometry, 
by interpreting the angle between the Wilson lines as rapidity. We discuss the relevance of our results 
for the phenomenology of near-forward hadronic processes at high energy, especially for processes 
with multiple color exchanges. 



I. INTRODUCTION 

Significant progress reached in the realm of non- 
perturbative QCD has been mostly related with ap- 
proaches based on the Euclidean formulation of the the- 
ory: numerical simulations using lattice gauge theory, in- 
stantons, monopoles, etc. By now, we know a great deal 
about the important or even dominant role of instanton- 
induced effects for correlation functions in a variety of 
hadronic channels, hadronic wave functions and form- 
factors, for a review see Unfortunately so far many 
of those results have not been translated to Minkowski 
space, a crucial step for understanding hadronic high- 
energy processes. It is however clear that there must 
be a very general and direct relationship between the 
hadronic substructure and the details of high energy reac- 
tions. Indeed, the non-perturbative modification of QCD 
vacuum fields induced by the valence quarks studied in 
Euclidean space-time should look like parton correlations 
in the transverse plane in a boosted frame. Many known 
features of partonic distributions, including spin and fla- 
vor of the sea quarks, point to their non-perturbative 
origin. Many more features (like fluctuations of these 
cross sections and correlations in the parton positions in 
the transverse plane which we briefly discuss at the end 
of the paper) are still to be studied in details. 

The first systematic step towards a semi-classical but 
non-perturbative formulation of high-energy scattering 
in QCD was suggested by Nachtmann [0, who has re- 
lated the scattering amplitude to expectations of pairs 
of Wilson lines. Semi-classical expressions with a similar 
pair of Wilson lines for DIS structure functions were also 
proposed by Muller ||: in contrast to their traditional 
interpretation as partonic densities, they were treated as 
cross sections for targets penetrated by small dipole-like 
probes at high energy. One systematic way to use these 
semi-classical expressions is to go back to the pertur- 
bative domain and try to improve on the diagrammatic 
approaches (like the celebrated BFKL Q re-summation) : 
see e.g. calculations of the anomalous dimension of the 
cross-singularity between two Wilson lines || or the anal- 
ysis of the path-exponents in Q . 



The approach we will follow in this paper is different: 
the Wilson lines in question are evaluated semi-classically 
using instantons. In order to be able to do so, one should 
start in Euclidean space-time, where those solutions are 
the saddle points of the functional integrals. The results 
are then analytically continued back to Minkowski space. 
Although it was not done before in this form, there are 
similar approaches in the perturbative context (e.g. JtJ 
and references therein). Another methodically close ap- 
proach to our analysis is ^|,[| where recent progress on 
the non-perturbative dynamics in N=4 SYM theory was 
used. In particular, the AdS/CFT correspondence have 
been used to evaluate the partonic cross section geometri- 
cally, using a deformed string in the curved anti-de-Sitter 
space. 

The instanton-induced processes to be considered in 
this work are either elastic scattering of partons, or 
quasi-elastic ones, with color transfer between them. 
They are very different from (and should not be con- 
fused with) multi-quanta production processes originally 
discussed in electroweak theory in connection with 
baryon number violation and later in QCD in connection 
to DIS [jll| . Such phenomena, associated with small-size 
instantons, are easier to evaluate and also they should 
lead to much more spectacular events. However, those 
lead to much smaller cross sections in comparison to the 
processes to be discussed below. 

In this paper, we will not aim at a development of a 
realistic model for high-energy hadronic reactions based 
on instanton physics. Instead, we will answer few ques- 
tions of principle, such as : Is it possible to assess non- 
perturbatively scattering amplitudes using the Euclidean 
formulation of the theory? How is the analytical con- 
tinuation enforced on the non-perturbative amplitudes? 
What is the magnitude of the instanton induced effects in 
comparison to the perturbative effects in the scattering 
of near-forward high energy partons? 

In section 2 we review the perturbative effects on 
the dipole-dipole potential, including the derivation of 
a renormalization group solution that can be tested us- 
ing QCD lattice simulations. In section 3 we extend the 
perturbative analysis in Euclidean space to the case of 
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scattering between two quarks and two dipoles. Partic- 
ular issues regarding the analytical continuation of the 
perturbative results to Minkowski space are discussed. 
In section 4, we discuss the effects of instantons on the 
static potentials for quarks and dipoles. At large dis- 
tances the results resemble perturbation theory apart 
from the large classical enhancement of (87r 2 /g 2 ) 2 « 10 2 , 
which is partially compensated by the diluteness fac- 
tor uqPq ~ (1/3) 4 |@ °f the instantons in the vac- 
uum. In section 5, we calculate the scattering amplitudes 
for quarks and dipoles in the one-instanton approxima- 
tion. The color preserving part of the amplitude is real 
and vanishes at high energy. The color exchange part 
is real but finite at high energy, thereby contributing 
to the near-forward inelastic scattering or re-scattering 
of partons. In section 6, we extend our discussion to 
two-instantons. We found that for two quarks the cross 
section is of the order of o~ qq ~ (noPo) 2 Pq, while for 
two dipoles it is further suppressed ~ "ggO^i^i/ Pa)- 
These results are supported by our calculations. In sec- 
tion 7, we discuss the possible role of instantons in cross- 
section fluctuations. Our conclusions and recommenda- 
tions are in section 8. 



II. PERTURBATIVE ANALYSIS OF 
POTENTIALS 

A. Dipole- dipole potential 

We start with the simplest analysis in Euclidean space, 
in which the perturbative expansion of two Euclidean 
Wilson lines leads to the well-known result for the po- 
tential between static charges. Indeed, by expanding two 
Wilson lines to first order in the gauge-coupling g, using 
the Euclidean propagator (A (x)A (y)) ~ l/(x — y) 2 with 
x, y located on two parallel but straight lines, and finally 
integrating over the relative time xo — yo, we readily ob- 
tain the Coulomb potential V(R) ~ a s (R)/R. 

Now, consider the case of the interaction between two 
color neutral objects, such as two static color dipoles. The 
simplest perturbative process in this case includes double 
photon/gluon exchange. The problem was solved in QED 
by Casimir and Polder who have shown that the 
potential at large distances R is 



where the polarizabilities a.1.2 are of the order of ai,2 ~ 
Tod 2 , d is the dipole size and To is some characteristic 
time (see below). This result differs from the Van-der- 
Waals potential 1/R 6 (valid at smaller R) because of the 
time delay effects. These observations were generalized 
to perturbative QCD in |l4j,[L5j. 

The Euclidean approach leads to the 7-th power of 
R in a simple way, provided that the following condi- 
tions are satisfied: (i) d\^ <C R which justifies the dipole 



approximation and identifies the relevant field operators 
(d-E) 2 ; (ii) both exchanged photons (or gluons) are emit- 
ted and absorbed at close xq and yo times. As a re- 
sult, the perturbative field correlator < E 2 (x)E 2 (y) >~ 
l/(x — y) 8 = l/(i? 2 + t 2 ) 4 , once integrated over the rel- 
ative time t = xq — yo, leads the result 1/R . 

The condition (ii) can be understood for complex sys- 
tems like atoms or hadrons in the following way: the 
first dipole emission excites the system from (usually an 
S-wave) ground state to (usually a P-wave) excited state, 
while the second dipole emission returns the system back. 
The energies of the intermediate state sets the character- 
istic life-time To ~ l/(Ep — E$). 

However for static dipoles the situation is different in 
QED and QCD. In QED the emission times of two ex- 
changed quanta are independent, but in QCD they are 
not. Even a static dipole can change its color degrees of 
freedom. Because different total color states of the dipole 
have different energies, thanks to the Coulomb interac- 
tion, we again have excited intermediate states. There- 
fore the characteristic time is determined by the differ- 
ence in Coulomb energy between the singlet and octet 
states 

1/t = AE = (3a s /2)/d . (2) 

Although the dipoles may be small d <C R, this time 
may still be long because in the perturbative domain the 
coupling constant is small g 2 (d) <C 1. As a result, there 
are two different regimes, when the distance R is large 
(i) R^S> To or small (ii) R <C To- In the former case again 
the power is 7 and the polarizability is[] a — 47rd 3 /3. The 
latter case is the Van-der- Walls domain. 



B. RGE analysis of the dipole-dipole potential 

On general grounds, the potential between two inter- 
acting dipoles can be shown to obey the following equa- 
tion 

where a s is the QCD running coupling and the averaging 
in (||) is carried in the presence of the two static dipoles 
a distance b apart. Generically, 

V(b) = V(b, a, fi, a.) w n (p a) K F{p 6, a.) (4) 

where p is the renormalization scale. Hence, 

£__!( ( . +1)V+ ,£) 



1 Amusingly, the result is just the volume of a sphere of radius 
d, from which the perturbative coupling constant g dropped 
out. 
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with (3 = da s /d\nfi is the QCD beta function. Inserting 
(§) into (|) yields 



where as usual s 
4m 2 and 



(pi + P2) 2 , t = (px- fci) 2 , s + t + u 



( K + l)V + b^- = A / d 3 x (TrF 2 {x)), , (6) W u (6) = P c exp (ig f + dr A(b 



+ Vl,2T) ■ Ul,2 



(14) 



which is the RGE equation satisfied by the dipole-dipolc 
potential. At large separations we may assume the 
dipole-dipole potential in quenched QCD to follow like 
a power law, i.e. 

V(6)«/i0io)»(M&F, ( ? ) 
turning (^) into an algebraic equation 

(1 + 7 + «) V(6) = A / d 3 x (TrF\x)) b . (8) 

Alternatively, the potential between two dipoles is a mea- 
sure of the energy density in the presence of two-dipolcs 



V(b) = / d 3 x (Qoo(x)) l 



(9) 



The combination of the RGE equation m) and the def- 
inition (^) yields a constraint between the exponents k 
and 7 in ([7]) asymptotically, namely 



7 = — 1 — K 



(3 1-R 
a s 1 + R 



with 



R 



Jd 3 x(B 2 (x)) t 
Jd 3 x(E 2 (x)) b 



(10) 



(11) 



a measure of the magnetic-to-electric ratio in the con- 
figuration composed of two static dipoles a distance b 
away from each other. For a self-dual field R = 1 and 
7 = — 1 — k if the asymptotic (|?|) is assumed. 



III. PERTURB ATIVE SCATTERING IN 
EUCLIDEAN GEOMETRY 

A. Quark-Quark scattering 

Generically, we will refer to quark-quark scattering as 

QA(Pi) + QB(p 2 )^Qc(ki) + Q D (k 2 ) (12) 

We denote by AB and CD respectively, the incoming 
and outgoing color and spin of the quarks (polarization 
for gluons). Using the eikonal approximation and LSZ 
reduction, the scattering amplitude T for quark-quark 
scattering reads |^,| 16 



AB.CD 



(s,t) w -2is / d 2 b e iq ^ b 



x((Wx(6) 



'AC 



(W 2 (0) 



I BD 



(13) 



The 2-dimensional integral in ( p^| ) is over the impact pa- 
rameter b with t = —q 2 ^, and the averaging is over the 
gauge configurations using the QCD action. The color 
bearing amplitude ( [13] ) allows for scattering into a sin- 
glet or an octet configuration, i.e. 



T = Ti 1 <g> 1 + T N 2^ (r a <g> r a ) 



(15) 



following the decomposition N c (£> N c = l(B (N 2 — 1). For 
gluon-gluon scattering the lines are doubled in color space 
(adjoint representation) and further gauge-invariant con- 
tractions are possible. For quark-quark scattering the 
singlet exchange in t-channel is + (pomeron) while for 
quark-antiquark it is - (odderon) as the two differ by 
charge conjugation. 

A quark with large momentum p travels on a straight 
line with 4- velocity x = v — p/m and v 2 = 1. In the 
the eikonal approximation an ordinary quark transmutes 
to a scalar quark. The argument applies to any charged 
particle in a background gluon field, with the following 
amendments: for anti-quarks the 4-velocity v is reversed 
in the Wilson line and for gluons the Wilson lines are 
in the adjoint representation. Quark-quark scattering 
can be also extended to quark-antiquark, gluon-gluon or 
scalar-scalar scattering. For quark-antiquark scattering 
the elastic amplitude dominates at large y/s since the 
annihilation part is down by \f—tfs. 

It can be described in Minkowski geometry in the 
CM frame with pi/m = (00SI17/2, sinh7/2, 0j_) and 
Pilm — (COSI17/2, — sinh7/2, OjJ with the rapidity 7 de- 
fined through cosli7/2 = y/s/2m. For s ^> m 2 the ra- 
pidity gap between the receding scatterers become large 
with 7 w log(s/m 2 ). The momentum transfer between 
the scatterers is q = p\ — k%, with g ~ 93 ~ t/\fs an d 
q 2 ^ = tu/(s — Am 2 ) rj —t. Hence q = (0,0, q±_) with 
q 2 = —q\ = t. Although the partons or dipoles change 
their velocities after scattering, this change is small for 
s ^> —t. This is the kinematical assumption behind the 
use of the eikonal approximation. 

In Euclidean geometry, the kinematics is fixed by not- 
ing that the Lorenz contraction factor translates to 



cosh 7 



v 1 



2m 2 



- 1 



(16) 



Scattering at high-energy in Minkowski geometry follows 
from scattering in Euclidean geometry by analytically 
continuing 9 — * —27 in the regime 7 w log (s/m 2 ) 3> 
1 . It is sufficient to analyze the scattering for p±/m = 
(1,0, 0_l), pil"m = (cos0,— sin8, OjJ, q = (0,0, q±) and 
b = (0,0, b±). The Minkowski scattering amplitude at 
high-energy can be altogether continued to Euclidean ge- 
ometry through 
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T AB ,CD(0 7 q) « 4m 2 sin 6> J d 2 b e iq ^ b 
x((W(9,b)-l) AC (W(0,0)-1) 



where 



W(6, 6») = P c exp (ig J dr A(b + vr) ■ v^j 



(17) 



(18) 



with v — p/m. The line integral in (|18j) is over a straight- 
line sloped at an angle 9 away from the vertical. 

In QCD perturbation theory, different time-ordering 
contributions to quark-quark scattering are shown in 
Fig. [j] to order g 2 . They contribute to the T-matrix as 
T = 27i + 2T 2 with (T -> oo) R 



71(0,6) 



cosfl 



(n - r 2 cos 9) 2 + r| sin 2 6> + b 2 



tantf 4tt 2 S \6 



(19) 



and 



W&) = ^ / dn 



o 

dr 2 

-T 

cos 9 



(n - r 2 cos 6>) 2 + rf sin 2 6> + 6 2 



(tt-0) g 2 / T 
tan0 4tt 2 ° S I 6 



(20) 



with 7? i 



) = — T\{tt— 9, b) as expected from geometry^. 
We note that the overall linear dependence in 9 reflects on 
the range of the gluon exchanged in rapidity space caused 
by our ordering in time. This dependence becomes 9 + 
(tt — 9) = it in the sum T, i.e. 



(21) 



as the ordering is unrestricted between and tt. All 
gluons between the spatial distance b and T are also 
exchanged, hence the infrared sensitivity of the quark- 
quark scattering amplitude in perturbation theory. This 



2 The color factors can be restored trivially. 

3 The reader may be puzzled by why we are emphasizing 
this simple point. We note that for more involved multi- 
gluon processes this cancellation is spoiled by color factors 
and powers of the angle survive in the answer: after the an- 
alytic continuation to Minkowski space these powers become 
powers of rapidity. They exponentiate and produce powers of 
the collision energy characteristic of Reggeon behavior (to be 
described elsewhere). 



sensitivity drops from the cross section (see below). We 
note that the order g 2 contribution to (|2l]) is of order 
s° after analytical continuation, in agreement with the 
general energy-spin assignment for vector exchange. We 
recall that the expected behavior is s J_1 for a spin-J ex- 
change. 

The contribution of ( pl| ) to T follows after integrating 
over the impact parameter b. The result in Euclidean 
geometry is 



T{9,q) =Am 2 sin 9 J d 2 be iq b T(9,b) 
g 2 Am 2 



-cos 9 '■ 



2 q 2 



dx Jq (x) logx 



(22) 



which can be translated into Minkowski geometry by an- 
alytical continuation through 9 — > —ij with q 2 = —t. 
In both geometries, T is purely real and divergent as 
t — > 0, leading to a differential cross section of the or- 
der of da/dt m g^/t 2 with a corresponding divergent 
Coulomb cross section a ~ g A /{—t m in)- In perturba- 
tion theory, the T matrix acquires absorptive parts and 
turns complex to higher-order, i.e. T = g 2 /t + ig i /t + ... . 
The Euclidean perturbative analysis can be carried out 
to higher orders as well, in close analogy with analytically 
continued Feynman diagrams . 





(a) (b) 

FIG. 1. One-gluon exchange between two receding par- 
tons, as discussed in Eq.(19) (a) and Eq.(20) (b). 



B. Dipole-Dipole Scattering 



We now consider dipole-dipole scattering 



D A (px) + D B (pa) -> D c (fci) + D D (fc 2 ) 



(23) 



emphasizing its color degrees of freedom. For simplicity 
we will assume both dipoles to have sizes d, and (in this 
section) average over their orientations. For pedagogical 
reasons, we start with a "naive" Euclidean approach at 
large impact parameter b, analogous to the calculation of 
the dipole-dipole potential above. This would be shown 
to lead to an incorrect answer for the high energy scat- 
tering amplitude. The reason will be given below along 
with the correct answer. 

We will assume that the impact parameter b is large 
in comparison to the typical time characteristic of the 
Coulomb interaction inside the dipole, i.e. !) > ro « 
d/g 2 . In the elastic dipole-dipole amplitude the dipoles 
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remain color-neutral, and may argue that the leading or- 
der is 2-gluon dominated. In analogy to the potential, 
one may rely on the Coulomb interaction inside the dipole 
to write the dipole-dipole effective vertex in the form 



S e ff - &e 



dr iy.iy F„ a F ua (x) 



(24) 



where the electric polarizability oe ~ (gd) 2 /£ with £ w 
g 2 /d its Rydberg energy |[4]]. (Higher order operators 
are suppressed by powers of the dipole size d.) 

In leading order in the dipole size, the scattering am- 
plitude then reduces to 

T(9, b) w as 2 xi fl ii u ±2\X2cr 

xf dndrJF^MF^F^)) , (25) 

J — oo 

with x\ = Win and xi = + The last expectation 
value can be unwound using free field theory to obtain 



T(0, 



1) al /ll 1 8 cos 2 6> 



7T° 



b e V 25 sin 5 sin 6 



(26) 



We note that the result ( pq ) diverges as — > 0. For the 
case = 0, we obtain the Casimir-Polder-type amplitude 



T(0,b) 



(A 2 - 1) Ta| 23 

7T 3 6 7 8 



(27) 



with T — * 00, which differs from the ^ by the occur- 
rence of the infrared sensitive factor T/b. 

The analytical continuation of (26) to Minkowski space 
shows that the first contribution is of order 1/s, while the 
second contribution is of order s. This implies that the 
total cross section is unbound, i.e. a ~ s, which is clearly 
incorrect. Indeed, on physical grounds the total cross sec- 
tion should be constant at large s. In Minkowski space 
it is easy to understand what went wrong. The electric 
field of a boosted dipole looks like a Lorenz contracted 
disk with a very small longitudinal width b/chy <C b. 
Clearly, at high energy the interaction time of two dipoles 
is of this order of magnitude, which is much shorter than 
the Coulomb time To. During this short time, the color 
rotation induced by the Coulomb interaction can be ig- 
nored. Therefore, the use of ( g4| ) in the form of a local 
2-gluon exchange is incorrect nT This point is actually 
missed in the Euclidean formulation as the Lorenz fac- 
tor is cos 9 1, Although any particular integral can 
be analytically continued from Euclidean to Minkowski 
space, kinematical approximations can only be inferred 
from the Minkowski domain where all parameters have 
their physical values. This will be understood through- 
out. 



So ignoring the Coulomb interaction and using the 
eikonal approximation, LSZ reduction and the analytical 
continuation discussed above, we can write the dipole- 
dipole scattering amplitude T in Euclidean geometry 
similarly to (O) with 



W(6,6) = i-Tr(P,<->q» 



ig l dr A(x) ■ v 

Co 



(28) 



where x is an element of Cg . In Euclidean geometry Cg 
is a closed rectangular loop of width d, that is slopped 
at an angle 9 with respect to the vertical direction. To 
leading order in the dipole-interaction, T can be assessed 
by expanding each Wilson-line ( pSj ) in powers of <?, and 
treating the resulting 2-gluon correlations perturbatively. 
The result is 



T(9,b) 



N 2 - 1 {gd) A cotan 2 



N 2 32tt 2 6 4 



(29) 



for two identical dipoles d\ = c?2 = d with polarizations 
along the impact parameter b. For small size dipoles, ( |29| ) 
is the dominant contribution to the scattering amplitude. 
The analytical continuation shows that cotan 9 — > — i, 
leading to a finite total cross section as expected. 



IV. INSTANTON EFFECTS ON THE 
POTENTIALS 

A. Generalities 

Instantons are self-dual solutions to the classical Yang- 
Mills equations in vacuum originally discovered in ref. 
p7[ . They are classical paths describing tunneling be- 
tween topologically inequivalent vacua of the gauge the- 
ory. In QCD, instantons were argued to be responsible for 
observable phenomena such as the resolution of the U(l) 
problem (large 77' mass ) fll8f a nd the spontaneous break- 
ing of chiral symmetry |!2]Jl9[ | . The interacting instanton 
liquid model (IILM) has been shown to reproduce multi- 
ple correlation functions, including hadronic spectra and 
coupling constants (for a review see [0). 

Instantons are also commonly used in other gauge the- 
ories, especially in supersymmetric gauge theories where 
supersymmetry makes their effects dominant in the non- 
perturbative regime. Indeed, some exact results (such as 
the effective low energy Lagrangian for N=2 supersym- 
metric theories derived by Seiberg and Witten and also 
the AdS/CFT correspondence suggested by Maldacena 
[] for the N=4 super-conformal theory) can be exactly 



4 Note that the result (^) is not based on this approxima- 
tion, and therefore is still valid. 



5 In fact, the 5- dimensional anti-de Sitter space emerges from 
the space of the instanton collective coordinates (the center 
position and size d 4 zdp/p 5 ) which will be extensively used for 
averaging below. 
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reproduced using exclusively the instanton calculus de- 
veloped in |po| . 

For the purpose of this paper the topology of instan- 
tons is not important: heavy quarks do not interact with 
fermionic zero modes, and high energy quarks for all 
purposes behave as heavy quarks. What is important 
instead is the following technical point: in the instan- 
ton field the path-ordered exponents can be evaluated 
analytically, since the color phase rotations take place 
around the same axis for a fixed path (the instanton is a 
hedgehog in color-space). The self-duality of the instan- 
ton field will also have an effect on some of our results. 
Once a path-ordered exponent is evaluated in the one- 
instanton field, the vacuum averages follow through the 
instanton ensemble average representing the QCD vac- 
uum (dilute phase). This includes averaging over the 
instanton center-position z M and size p. Specifically, we 
will use the measure 



dn = dpd z 



4 D(p) 



(30) 



for both instantons and anti-instantons. The integral 
over z can be sometimes carried out analytically, but 
most of the times will be done numerically. The un- 
derstanding of the instanton size distribution D(p) re- 
mains an open problem. Naive semi-classical results sug- 
gest@: 



D (p) « C 



N, 



(3T77t) ex P(" 



9 2 (p) 
(pA)( u /3)^ 



8ir' 
9 2 (p) 



-(2/3)N f 



(31) 



where CV c is a constant depending on the number of 
colors N c . We have used the asymptotic freedom for- 
mula in the exponent to show that this density dramat- 
ically grows with the instanton size p. However in the 
true QCD vacuum instantons and antiinstantons inter- 
act with each other and othe quantum fields, so that the 
real function D(p) deviates from the semiclassical one for 
large sizes. 

For qualitative estimates we will often use parameters 
of the instanton liquid model [O] , which assumes that all 
instantons have the same size 



dn(p) — nod zdpS(p — po) 



(32) 



where no is the total instanton (plus anti-instanton) with 
a typical radius po, i.e. 



n 



1 fm A ; po » 1/3 fm 



(33) 



These values were deduced from phenomenological data 
extracted from the QCD sum rules, the topological suc- 
ceptibility and the chiral condensate long before direct 
lattice data became available. In Fig. we show a sample 
of such lattice measurements, together with the parame- 
terization for the instanton suppression suggested in pl| . 
Specifically, 




FIG. 2. (a) The instanton density dn/dpd 4 z, [fm -5 ] versus 
its size p [fm]. (b) The combination p~ 6 dn/dpd 4 z, in which 
the main one-loop behavior drops out for N c — 3, Nf — 0. 
The points are from the lattice work [Q , for this theory, with 
/3=5.85 (diamonds), 6.0 (squares) and 6.1 (circles). Their 
comparison should demonstrate that results are rather lat- 
tice-independent. The line corresponds to J21|. 



dn(p) = dno{p)e 



-2irap 



(34) 



can be used for averaging in any integral over the in- 
stanton density. Typically, the string tension a « 
(0.440 MeV) 2 , so that 



< P 5 >r 



(0.28 fm) 2 ; 
(0.31 /m) 4 ; 
(0.32 fm) 5 , 



(35) 



which shows that the difference between the realistic av- 
erages and simple powers of po are relatively small. We 
will ignore these differences below. 

In the analysis to follow, the parameters capturing the 
instanton physics will appear as two dimcnsionless quan- 
tities: (i) a small diluteness parameter and (ii) a large 
action of an instanton (per h): 



So 



9 2 (po) 



(10-15) (36) 



The small factor is a penalty for finding the instanton, 
and the large factor is a classical enhancement relative 
to perturbation theory. Their interplay would cause par- 
ticular effects to be parametrically large or small. 



B. Static quarks 

At the one instanton level, the various potentials for a 
static quark-antiquark potential have been assessed long 
ago ^3|, including the spin-dependent part. We will 
briefly review this assessment for completeness. We re- 
call that the various components of the potential follow 
from the rectangular T x R Wilson loop 



G 



V(R) = -- lim In < W(T,R) > 

1 T — >oo 



(37) 



evaluated in a classical instanton field, after averaging 
over the instanton position. In the Wilson loop, the 
path-ordered exponents Pexp(ig J A^dx^) can be eval- 
uated analytically as the instanton locks the color ori- 
entation to space. Indeed, the static potentials involve 
Aq ~ Vo u( x ~ z )v ~ ( x — z ) a where (x — z) refers to the 
distance between the quark position and the 3-d coordi- 
nate of the instanton center^. The resulting color rota- 
tion angle a p3| and the unit vector n a around which 
the rotation takes place are defined through 



W = cxp — in 



r a {z a -r a ) 



((r a - 
exp(—iTTT a n a a) 



p 2)l/2 



(38) 



If all relevant distances are comparable, \r a — z a \ ~ po, 
the rotation angle is 0(1), showing that the expansion 
in field strength is in general not justified. For a small- 
size dipole, the potential is small V(R — + 0) w R 2 , since 
the path-ordered lines in W are close enough to cause 
partial cancellation. However when R w p Q , and both 
path-ordered lines happen to be on the opposite sides of 
the instanton center, the color rotations on both lines 
adds up and the potential becomes roughly linear in R 
and more sizable. Finally, when the dipole is too large, 
the potential saturates^ 

The quark-antiquark potential calculated in ^3| can 
be expressed as 



(39) 



V(R) = J dn(p)p 3 F(R/p) 

where the dimensionless function F is defined as 

^ = /0Tr(l-W lW t). 

The trace-part of the integrand is 

2(1 — cos ai cos 02 — n% ■ ni sin ct\ sin ol-i) . 



(40) 



(41) 



where the angles ai and vectors nf are defined in (38). 
This function is shown in Fig. ^| a. In order to empha- 
size the small-R "dipole limit" V(R) ~ R 2 (to be im- 
portant for what follows), we have also plotted the ratio 
of this function to its dipole limit in Fig. || b. One can 
see that the dipole approximation has an unexpectedly 
large range of applicability: this ratio does not change 
appreciably (less than 25%) till Rm po. One may expect 



The time position of the instanton 24 is irrelevant. 
7 In ref. j3^] one of us has noticed that this behavior is surpris- 
ingly similar to that experimentally observed in deep inelastic 
scattering, if Q 2 dependence of structure functions is treated 
as dependence of the cross section on the dipole size. 



similar accuracy of the dipole approximation in other ap- 
plications to be discussed. 




FIG. 3. (a) F as defined in (J3S|) as a function of the 
quark-antiquark distance R measured in units of the instanton 
size po; (b) the rescaled function F(R/ pa)/ R 2 to exhibit the 
accuracy of the dipole approximation at small R. 

For large R the potential [0 goes to a constant plus 
a Coulomb term 



V { R^oo) = v(%D{ p) -% d £ D{p ) 
J p l 3R p 



(42) 



which can be interpreted as the instanton contribution 
to the mass and charge renormalization, respectively. It 
is instructive to compare the magnitude of the latter to 
the perturbative potential, through 



Vvnst 7T , 4 87T 2 
Vpert * 9 (PO) 



(43) 
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with Vp er t = Aa s /2>R. The ratio is the product of the di- 
luteness parameterf\ (the fraction of space-time occupied 
by instantons) times the classical enhancement through 
the instanton action (per K). Using the phenomenologi- 
cal parameters discussed above, we observe that the di- 
luteness is compensated by the classical enhancement, so 
that the instanton corrections at R ~ po are actually 
comparable to the perturbative Coulomb effect. 

However, instantons are not the only non-perturbative 
effects contributing to the static quark-antiquark poten- 
tial. At large R confinement in the form of a QCD 
string with V con f ~ oR dominates. In fact, already 
for R ps po ~ 0.3 fm confinement is dominant, with 
the instanton-induced potential accounting for only 10- 
15% H For a detailed study of these issues at the multi- 
instanton level, one can consult refs pl| for a numerical 
analysis and [E5| for analytical results. 




C. Static dipoles 



Unlike the quark-antiquark potential, the dipolc-dipole 
potential is insensitive to confinement, and the instanton- 
induced interaction may be easier to identify. In the lat- 
ter case, we will consider two cases where the characteris- 
tic time within the dipole is either (i) short tq ~ d/g 2 -C 
p or (ii) long tq 3> po in comparison to the instanton 
size. These two cases translate to a magnitude of the 
dipole field Aq ~ g/d which is large (i) or small (ii) in 
comparison to that of the instanton field ~ 1 j gp. 

In the case (i) the static potential can be written in 
terms of the polarizabilities, and the correlator of gluo- 
electric fields 



V(R) = a ia2 / dr < E 2 (t, R)E 2 (0, 0) > 



(44) 



This field strength correlator can be evaluated by substi- 
tuting the expression for the instanton field 



E 2 {x) =B 2 {x) 



96p 4 



1 



g 2 {{x-z) 2 + p 2 Y 



(45) 



The averaging of the correlator over the location of the 
instanton position z can carried out analytically |27[ 



< (.9G^(x)) 2 ( 5 G^(0)) 2 >= 



(46) 



where the last term was added to account for the pertur- 
bative contribution. The dimensionless function describ- 
ing the instanton contribution is 



8 The coefficient in front of 7r 2 p 4 /2 happens to be the volume 
of a 4-sphere. 

9 The claim made in [Q, that instanton effects account for 
the confining potential is incorrect. 




FIG. 4. (a) Field strength correlator Ui ns t as defined in 
(0) versus the distance in units of the instanton size x/po- 
(b) Ratio of the instanton-induced term in the correlator to 
the perturbative one, versus the distance. 



Uinst (y) 



12288tt 2 



4)5 



y 6 (y 2 

x(y 8 + 28y 6 
240 



94 y 4 - 160 y 2 - 120 



y 



x (y 6 + 2 y 4 + 3 y 2 + 2) arcsinh(y/2)) . (47) 

Its behavior is shown in FigJJa. Its ratio to the pertur- 
bative contribution to the same correlator (for g = 2 
or a s = 0.32) is shown in fig.|]b. As expected, it is 
small at small distances i < po. At large distances, 
the instanton-induced contribution has the same behav- 
ior Hi ns t ~ 1/-R 8 j as the perturbative one. Furthermore, 
the ratio of the two is about 30, much more than the 
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"instanton-induced charge renormalization" (^3|) we dis- 
cussed in the preceding subsection. About the same is 
found in the potentials themselves (the correlator inte- 
grated over the time difference) as shown in Fig. || d. 
The perturbative behavior is dominated by two gluons 
rather than one, and therefore the instanton effect oc- 
curs with a classical enhancement squared: 



Vi, 



v, 



pert 



(noPo) 



9 2 (po) 



(48) 



This feature implies that instanton effects are much more 
important for dipole-dipole interactions at R « po ~ 0.3 
fm than the perturbative Casimir-Polder effects. We will 
argue below that this is generic for all processes demand- 
ing multi-gluon exchanges, and that instanton-induced 
processes can become dominant in this case. 

In the case ii), the dipoles can be considered quasi- 
static in time, To ~ d/g 2 3> po, and the time evolution 
of the color degrees of freedom due to the Coulomb in- 
teraction can be ignored. In other words, the dynamics 
is driven entirely by the instanton field. The potential 
between two dipoles is now 



V dd {R) = / dn(p)p 3 F dd {R/p) 



with 



F, 



dd 



d 3 Z 

NoP : 



■(l-TrWiTrW 2 ) 



(49) 



(50) 



Here W arc rectangular Wilson loops for each dipole, 
traced separately. Averaging over the instanton position 
can be done numerically. The results are shown in Fig|3. 
The outcome is proportional to d\d\ (dipole-moments) 
rather than ot\ oli (electric polarizabilities), when d is 
reasonably small in comparison to po. The large distance 
potential is few % that of V(R) « dfdlpl/R 7 . Note that 
it is larger than the perturbative one since p\ is assumed 
to be much larger than g?ig?2, but both answers have the 
same (zeroth) power of g. 

In general the dipole-dipole potential cannot be ap- 
proximated by the correlator of scalars E 2 , as can be 
checked through its dependence on the relative orienta- 
tion of the dipoles. Even in the dipole (quadratic) ap- 
proximations for sufficiently small dipoles (di <C p) one 
can define 4 invariant functions for the dipole-dipole in- 
teraction 

V(R) = d\d\d\d™ {A(R)S ij S lm + 

\B{R){n l n^ lm + n l n m 8 l3 ) + 



C(R)n i n l 8 jm + D(R)n l r^n l n T 



(51) 



The first function A(R) accounts for the spin-zero gluonic 
operator E 2 discussed at the beginning of this subsection. 
However, as one can see from fig(|j|), other functions also 
contribute. In (a) we compare the xx orientation (or 



A+B+C+D) with the xy one (or A+B/2) and see a clear 
difference. In (b) we note the dependence on the rotation 
angle for one of the dipoles, which shows a clear cos 2 9 
behavior expected from the expression above. 
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FIG. 5. (a)Two correlated Wilson loops as a function of 
the distance R between their centers, divided by d\d% for two 
dipole sizes, di = di = 0.1 (circles) and 0.4 po (squares). The 
agreement between the points show that the dipole scaling 
holds well for such sizes. Also two dipole orientations are 
shown. The open points are for both dipoles oriented in the 
x direction (the same direction as the distance between them 
R) while the closed points for the xy orientation. The dis- 
agreement between those means that the dipole-dipole forces 
depend on the orientations. Further details on the orientation 
dependence are shown in Fig.(b). Fdd(R = po) shows the de- 
pendence on the orientation angle of one dipole in the x-y 
plane. The solid points and curves are for a 4-d lattice-type 
integration over the instanton center, and the open points are 
for the alternative Monte-Carlo integration: their spread from 
the curve shows the magnitude of the uncertainties involved. 



V. ONE- INSTANTON EFFECT ON SCATTERING 

A. Quark-Quark scattering 

Our first step now is the generalization of ( |38|) to an 
arbitrary orientation 8 of the Wilson line. The analytical 







continuation to Minkowski space follows from 9 — > iy 
with y identified as the rapidity difference between the 
receding partons. The untraced and tilted Wilson line in 
the one-instanton background reads 



where 



W(0, b) = cos a — it ■ h sin a 



b „b 



(52) 
(53) 



the high energy limit) changes sign, before decreasing as 
a power law to zero at large b. 



and a — irj/ 'y/'y 2 + p 1 with 



n • n 



(z4sin# - z 3 cos#) 2 + (b - zj_) 2 



(54) 



The one-instanton contribution to the untraced QQ- 
scattering amplitude follows from the following correlator 



(W AC (9,b)W BD (0,0)) 



n 



d 4 2 



x (cosa cosa Iac Ibd 

- N 2_ 1 " ' nsinasina (t")a C (t") bd ) 



(55) 



where the (under) bar notation means the same as the 
corresponding un-bar one with 9 = and 6 = 0. 
Furthermore, 



^Tr(W(^,6)W(0,0)) 



2n f 4 ... , 

-jy- J a z (cos a cos a — n • n sin asm a; 



(56) 



The integrand in ( |56| ) can be simplified by changing vari- 
able (Z4 sin^— Z3 cos 9) — > z± and dropping the terms that 
vanish under the z-integration. Hence 



±Tr(W(M)W(0,0))^ = ^ J 



V sin0 



cos a cos a 



tan ( 



■ sin a sin a 



zj_ ■ b 



77 



(57) 



The tilde parameters follow from the un-tilde ones by set- 
ting 6 = 7r/2. We note that 7 = 7= \z\. After analytical 
continuation, the first term produces the elastic ampli- 
tude which decays as 1/s with the energy. The second 
term corresponds to the color-changing amplitude, ft is 
of order s° and dominates at high energy. Specifically 



^Tr(W(tf,fo)W(0,0))) = 



2no f 1 
N c Vsinfl 



F cc {b/p ) 



tan 9 



F ss (b/p ] 



(58) 



We show in Fig. p| t he numerical behavior of the two 
contributions in (|58|). Note that the second function 
(which describes color-inelastic collisions and survives in 



FIG. 6. (a) and (b) show the two functions noF cc and 
noFss defined in (^), versus the impact parameter b (in units 
of the instanton radius). 



B. Dipole-Dipole and Multi-parton Scattering 

One can directly generalize the calculation of the 
quark-quark scattering amplitude to that of any number 
of partons. For that, we assume that they all move with 
high energy in some reference frame and opposite direc- 
tion: in Euclidean space those would propagate along two 
directions, with parton numbers N% and N2 respectively. 
Any one of them, passing through the instanton field, is 
rotated in color space by a different angle a, around a 
different axis Hi, depending on the shortest distance be- 
tween its path and the instanton center. Integration over 
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all possible color orientations of the instanton leads then 
to global color conservation. 

Before discussing specific cases in details, let us make 
a general qualitative statement about such processes. 
We have found in the previous section that (the color- 
changing) quark-quark instanton-induced scattering has 
a finite high energy limit. For perturbative n-gluon ex- 
change a factor of a™ is paid, while for an instanton me- 
diated scattering a factor of noPo 1S P a id ( tne P r i ce to find 
the instanton at the right place), no matter how many 
partons participate. Since the instanton vacuum is dilute, 
the one-gluon mediated process dominates the instanton 
one. However, the situation dramatically changes for two 
or more gluon exchanges: the instanton-induced ampli- 
tude is about the same for any number of partons, pro- 
vided that all of them pass at a distance w po from the 
instanton center. 

Now, consider a dipole configuration of size d chosen 
in the transverse plane of a qq located on a straight-line 
sloped at an angle 9 in Euclidean space. Let AA be the 
initial color of the dipole and CD its final color. The 
Wilson loop with open color for the dipole configuration 
in the one-instanton background is 

W^a (0> b) — cos a_ cos a + 1 C d 
-Mcos a_ sin a + R afc n+ (T a ) DC 
-zsina_ cosa + R afc h b _ (r a ) DC 
+sina_ sina+ R ah R cd h b _ {T c T a ) DC ■ 

We have defined 

7T7± 



dipole-dipole scattering amplitude through a single in- 
stanton is 



(59) 



a± 



V7± + P 2 



7 ± = (z 4 sin# — z 3 cos9) 2 + (z± 



b± - 2 

2> 



n = (z 4 sin<9 - z 3 cos0) 2 + (b - z ± ) 2 - — (60) 



with n± • n± = j±. The scattering amplitude of an ini- 
tial dipole through an instanton after averaging over the 
global color orientation R is 

2 

— (cosa_ cosa+ + n_ ■ n+ sino;_ sina+) Icd (61) 

^ * c 

which reduces to the color-singlet channel. Specifically, 



W(6, b) — — (cos a_ cos a+ + n_ ■ n + sin a_ sin a+) 



The dependence in (^31-62) can be readily eliminated 
by carrying the integration over the instanton position 
z through the same change of variable discussed in the 
quark-quark scattering, resulting in an amplitude that 
depends only on l/sin#. In Minkowski space this trans- 
lates to 1/s which vanishes at high energy. Indeed, the 



(W(0, 6) W(0, 0)) 



rip 
sin 9 



d 4 zW(9,b)W{0,0) (63) 



where VV follows from W by setting 9 = ir/2. Note that 
in this case W(0, 0) = VV(0, 0). 

It is clear from ( p9| ) that while scattering through an 
instanton, the dipole has to flip-color to keep track of the 
velocity of the quarks in the dipole. The process is color- 
inelastic and therefore only contributes to the inelastic 
amplitude to first order in the instanton density no, and 
to the elastic amplitude to second order in the instanton 
density, a situation reminiscent of one- and two-gluon 
exchange. 

The dipole-dipole scattering amplitude with open-color 
in the final state can be constructed by using two dipole 
configurations as given by ([39]) with a relative angle 9. 
After averaging over the instanton color-orientations we 
obtain 

W^(0,&)W##(O,O) = 

J- Wi 1 C D 1cd> + 1 Wjv c 2 -i (r a ) D c {T a ) D >c> , 



(64) 



with the singlet part 

Wi = cosa_ cosa+cosa_ cosa + 
+n_ • n + n_ • n + sino;_ sina-f sina_ 
+n_ • n + cos a_ cos a+sin a_ sin a + 
+n_ • n + sina_ sina+cosa_ cosa + 

and the octet part 



(65) 



Wjv2_! = — cosa_ sina+cosa_ sina + n + ■ n + 
— sina_ cosa+sina_ cosa + n_ • n_ 
+cosa_ sina+sina_ cosa + n+ • n_ 
+sino!_ cosa+cosa_ sina + n_ • n + 
— sina_ sina+sina_ sina + 

x(n_ • n n + • n + n • n + n + • n ) . (66) 

As in the case of quark-quark scattering, the (color) elas- 
tic dipole-dipole amplitude scales as 1 /sin 9 and vanishes 
at high energy after analytical continuation. However 
the (color) inelastic part of the amplitude is not. After 
performing the shift of variables described before, the 9 
(62)lependence drops from all the angles a. There is a re- 
maining 9 dependence in the four combinations n • n. In 
general, the 9 dependence in the latter is linear in s'm9 
or cos 9, and one may worry that the last term in (|6^ ) 
may involve higher powers of the trigonometric functions, 
which would yield to an unphysical cross section growing 
as s after analytical continuation. We have checked that 
this is not the case, since 
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n + n + • n — > d (z\ — cos 9 z^z'^ 



where z\ is the new Z4 after the change of variable. More- 
over, the cos# term drops in the integral over z (odd 
under Z3 — > —23), making this contribution to (|66| ) sub- 
leading at high-energy after analytical continuation [^]. 
Finally, we note that all sin6* contributions in ([36]) drop 
following similar parity considerations. As a result, the 
pertinent octet contribution to the scattering amplitude 
is proportional to cotan0 which is l/itany = 1/iv after 
analytical continuation. 

We have assessed numerically the function 



F N 2_ 



b d 
Po' Po 



d 4 zW N 2_ 



(67) 



which is shown in Fig. for different dipole sizes. We find 
that the dipole approximation scaling F N 2_ 1 ~ d 2 works 
well, even for sizes as large as the instanton size d = pQ. 
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FIG. 7. (a) Fjva.iXM) defined in @ versus the 
dipole-dipole impact parameter b (in units of the instanton 
size po). Each curve corresponds to a different dipole size d 
(same units). 



VI. TWO- INSTANTON EFFECT 

We have shown above that the instanton contribu- 
tion at large s but small t behaves in a way similar to 
one-gluon exchange: only color-inelastic channels sur- 
vive in the high energy limit. This means that the 
contribution to the total cross section appears in the 



amplitude squared, leading naturally to the concept of 
two-instanton exchange. The latter contribution to each 
Wilson-line is more involved. To streamline the discus- 
sion we will present the analysis of the two instanton con- 
tribution to the differential cross-section of quark-quark 
scattering at high energy. Similar considerations apply to 
dipole-dipole scattering as we briefly mention at the end 
of this section. Indeed, for the quark-quark scattering, 
unitarity implies that the two-instanton contribution to 
the differential cross section is 



da 
lit 



BD 
AC 



(68) 



CD 



with the averaging over the initial colors A, B under- 
stood. Inserting ( |l3| ) after the substitution we ob- 
tain 



da 
~dt 



with 



4no 



dbdb' e lq {b ~ b,) I J 



1 



(N 2 - 1) 



K 



(69) 



J = y d 4 2 (cosa — l)(cosa — 1) 

x I d 4 z'(cosa' — l)(cosa' - 1) 
K = / d 4 z n • n sin a sin a 

a z n • n sin a sin a . 



(70) 



The primed variables follow from the unprimed ones 
through the substitution b — > z',b'. For large ^/s, 
J~(l - F cc )(l ~ F' cc )/s 2 Q and K = F ss FL, so that 



16 n\ 



da 



dt N 2 C {N 2 ~1) 



dbe lqh F s 



ss J- ss ' 
2 

Po 



(71) 



In particular, the forward scattering amplitude in the 
two-instanton approximation is 



a(t = 0) 



16 nl 



N 2 (N 2 - 1 
dq\ 



dbe^F^A — 



. (72) 



which is finite at large \fs. Hence, for forward scattering 
partons in the instanton vacuum model, we have 



a qq w (n po) 2 pi 



(73) 



10 This cancellation is not generic. Indeed, the square of this 
contribution would be leading. 



Up to self-energies. 
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Clearly, the present analysis generalizes to the dipole- 
dipole scattering amplitude by using (^J) instead of (|5^ 
and proceeding as before. The outcome is a finite scat- 
tering cross section, 



a(t = 0) 



4n 2 



dqj 



dbe tq - h F l 



A - 



b d 

Po' Po 



(74) 



Generically, the dipole-dipole cross section relates to 
the quark-quark cross-section in the forward direction 
through 



0~dd 



(did 2 



pi 



(75) 



It is instructive to compare our instanton results to 
those developed by Dosch and collaborators pa] in the 
context of the stochastic vacuum model (SVMjTln brief, 
in the SVM model the Wilson-lines are expanded in pow- 
ers of the field-strength using a non-Abelian form of 
Stokes theorem in the Gaussian approximation. A typi- 
cal hadronic cross section in the SVM model is 



< {gGf > 2 



,10 



F (R h /a) 



(76) 



where the first factor is the "gluon condensate" , a is a 
fitted correlation length, F is some dimensionless func- 
tion depending on the hadronic radius R^. Although our 
assumptions and those of f28f| are very different regard- 
ing the character of the vacuum state, it is amusing to 
note the agreement between ([75]) and (|76|). Indeed, the 
correlation length a of the SVM model is related (and in 
fact numerically close) to our instanton radius po ~ 1/3 
fm, while the gluon condensate < (gG) 2 > of the SVM 
model is simply proportional to the instanton density no 
in the instanton model. 

The most significant difference between these two ap- 
proaches apart from their dynamical content and the way 
we have carried the analytical continuation, is the fact 
that we do not expand in field strength. In fact, in the 
instanton model there is no parameter which would al- 
low to do so for strong instanton fields. This difference is 
rather important as it is on it that our conclusion regard- 
ing multiple color exchanges is based. (In the SVM model 
with Wick-theorem-like decomposition, those would be 
just products of single exchanges, like in pQCD.) 



VII. CROSS SECTION FLUCTUATIONS 



quantum mechanical states [^. Hence, the quantum sys- 
tem is characterized by some amplitude of probability 
through its wave function, and its corresponding scatter- 
ing cross section is probabilistic with a probability distri- 
bution -P(cr). This idea was originally suggested by Good 
and Walker |29(| , who emphasized that inelastic diffrac- 
tion is a way to quantify this distribution via the second 
moment Act 2 =< (ct 2 — < a >) 2 >. 

The extraction of this and the next (cubic) moment for 
the pion and the nucleon using available data has been 
carried out years later [ p0| allowing for a reconstruc- 
tion of the distribution P(a)- A striking aspect of these 
results is that the nucleon fluctuations are large and com- 
parable to the pion fluctuations. This outcome does not 
fit with the constituent quark model where the pion is a 
2-body system, and the nucleon is a 3-body system, with 
more degrees of freedom. One of us [ [jlf had already no- 
ticed that this can be a further indication for strongly 
correlated scalar diquarks in a nucleon. An experimen- 
tal test for this idea is to measure cross section fluctu- 
ations for a decuplet baryon such as In the latter 
there are no diquarks, and smaller fluctuations (typical 
of a 3-body state) are expected. Another aspect of these 
fluctuations worth mentioning here is that they seem to 
be maximal for y/s ps 100 GeV, decreasing at very large 
energies. It supports well the idea that the "most fluctu- 
ating" partons are at x ~ 10~ 2 , while at much smaller x 
one basically approaches a non-fluctuating black disk. 

Although in the present paper we have limited our dis- 
cussion to issues of methodology, it is worth pointing out 
that the present concept of fluctuations in cross section 
can be used to discriminate between the instanton effects 
herein described and other descriptions based either on 
perturbative multi-gluon exchange or non-perturbative 
vacuum structures. 

Indeed, the standard multi-photon exchange in QED 
leads to an (eikonalized) exponential scattering ampli- 
tude, with Poisson-like fluctuations. If the mean-number 
of quanta exchanged < n >^> 1 (e.g. for heavy ions with 
large Z« 1/a), the distribution becomes narrow and we 
approach a classical limit, with weakly fluctuating scat- 
tering. Modulo color factors, the same conclusion ap- 
plies to multi-gluon exchange in QCD. In contrast, the 
instanton-induced effects have completely different statis- 
tical properties. The field of the instanton itself is clas- 
sical, hence coherent. However, the distribution over the 
instanton size and position is quantum (in contrast to the 
Coulomb field of the ion just mentioned), thereby lead- 
ing to cross section fluctuations. The latter arc further 
enhanced by the diluteness of the instanton ensemble: 
the quark may appear very black, provided a tunneling 



In so far, we have considered the average value of the 
cross section for a parton in a state of unit probability. 
However, partons and in general hadrons, are complex 



12 A truly elementary particle may have only one state and 
non-fluctuating cross section: it may have diffraction but no 
inelastic diffraction. 
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event happens to be close to it, and rather transparent 
otherwise. As noticed already in fi~2]| , quarks are "twin- 
kling" objects, as the associated gauge/quark fields are 
strongly fluctuating. 

To quantify some of these statements we show in 
Fig.(||) how such distribution looks like. We plot \F ss (b — 
1)| 2 , at fixed impact parameter b=po. The distribu- 
tion corresponds to instantons being homogeneously dis- 
tributed in the 4d sphere around the center of the colli- 
sion point, with a radius R s « 2.2p such that n 2 R^/2 = 
1/riQ, or in a smaller sphere within R < p . However 
the resulting amplitude is highly inhomogeneous, with a 
large peak at small amplitude, and a long tail at large 
amplitude. Comparing the solid and dashed curves, one 
can see that the latter is due to instantons sitting near 
the center of the system. 




0.2 0.4 0.6 0.8 1.2 

Fss(b=1) 2 



FIG. 8. The distribution \F ss (b = po)\ 2 with instantons 
filling homogeneously the 4d 'Wigner-Seitz sphere' of radius 
2.2 po (solid histogram) or the smaller sphere of radius 1.0 po 
(dashed histogram). 

VIII. CONCLUSIONS AND OUTLOOK 
A. Conclusions 

Several new instanton-generated phenomena have been 
studied in this work: static potentials for color dipoles, 
and high energy quark-quark and dipole-dipole scatter- 
ing. The nature of the instanton effects makes their con- 
tribution to these processes different from the contribu- 
tion expected in perturbation theory. 

Overall, the magnitude of the instanton contribution 
is governed by two competing factors: (i) a diluteness 
factor uqPq <C 1 reflecting the fact that their density in 
the QCD vacuum is small (noPo "C 1), and (ii) a clas- 
sical enhancement factor, the instanton action which is 
large (So/% s» 10 ^> 1). Naturally, the more partons 



are involved in a particular process, the more powers 
of a s appear in the perturbative result for a particular 
process. This penalty does not apply to the instanton 
contribution. One way to quantify this difference is to 
note that the ratio of the instanton-to-perturbative con- 
tributions contains a power of the classical enhancement 
parameter, and this power grows with the number of par- 
tons involved. Typically the first power due to the clas- 
sical instanton enhancement cannot really compensate 
for the small diluteness of the instantons in the vacuum. 
However, the second power is already sufficient to make 
the instanton effects larger than the perturbative ones as 
we have now established for the potentials. Indeed, the 
dipole-dipole instanton-induced potential exceeds signif- 
icantly (by a factor ~ 25) the perturbative contribution 
for distances R > p®. 

Based on these ideas, we have extended the analysis to 
near-forward parton-parton scattering amplitudes, treat- 
ing in details the case of quark-quark and dipole-dipole 
scattering. Key to our analysis was the concept of analyt- 
ical continuation in the rapidity variable, which we have 
applied to both the perturbative and instanton analysis 
for comparison. 

In the perturbative analysis, one- and two-gluon ex- 
changes differ fundamentally in the sense that the for- 
mer is color-changing (inelastic), while the latter is color- 
preserving (elastic). Indeed, the two-gluon exchange 
mechanism |32|| constitutes the starting ground for the 
soft pomeron approach to dipole-dipole scattering. Since 
the instantons can be viewed as multi-gluon configu- 
rations (classical fields), we have suggested that they 
maybe a viable starting point to analyze soft parton- 
parton scatterings. We have shown that the instanton- 
induced amplitudes involve also color-elastic and color- 
inelastic channels. After analytical continuation, the 
one-instanton contribution to the color-elastic channel is 
purely real and vanishes as 1/ yfs (much like a scalar ex- 
change). In other words, in this work a single instanton 
is not "cut" , its multi-gluon content is not used. Instan- 
tons contribute to soft parton-parton scattering like the 
t-channel gluons mostly through color exchange chan- 
nels, or through re-scattering in the elastic channel. The 
leading instanton contribution involves a two-instanton- 
prong channel, and yields a finite elastic parton-parton 
scattering amplitude after analytical continuation in ra- 
pidity space. Our result is reminiscent of the one reached 
in the stochastic vacuum model [ pi) , although our as- 
sumptions and methodology are different. 

B. Outlook 

The results we have derived were achieved in Eu- 
clidean space prior to our pertinent analytical continua- 
tion. Therefore, they are testable from first principles by 
repeating our analysis using instead lattice QCD simula- 
tions. Indeed, the non-perturbative dipole-dipole forces 
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could be studied. In contrast to the quark- antiquark po- 
tential and to the best of our knowledge, those forces 
have not been investigated on the lattice. Also, the var- 
ious scattering amplitudes discussed in the present work 
can and should be looked at, leading to multi-parton am- 
plitudes as we have qualitatively discussed. Note, that 
not only the potentials and scattering amplitudes them- 
selves can be derived, but the degree of their correlation 
with the presence of instantons in the underlying config- 
urations can be revealed as well, using lattice techniques 
such as "cooling" and alike to help discriminate instan- 
tons by their topological charge. 

Regarding the applications of our results, we admit 
that there remains a significant distance to the descrip- 
tion of real hadronic processes. Although we hope to 
cover further phenomenological applications elsewhere, 
we still would like to comment on two broad but impor- 
tant dynamical issues: (i) the mechanism of color rear- 
rangements in high energy collisions and (ii) the issue of 
hadronic substructure of the non-perturbative effects in 
the hadronic wave functions. 

It is generally accepted that high energy hadronic pro- 
cesses can be split into three stages: (i) formation of 
hadronic wave function (to which we turn later); (ii) color 
re-arrangements of partons in a collision; and (iii) de- 
cay of the arising system into multi-hadron final states. 
It is further believed that at stage (iii) color flux tubes 
(QCD strings) are formed with basically unit probabil- 
ity^, so that one can ignore them in the calculation of 
the cross section. Such assumption is implied in any 
perturbative approach (such as the Low-Nussinov gluon- 
exchange model |p2[ ), and we assume the same is true for 
instanton-induced color exchanges as well. 



1 The fluxes are described by multiple phenomenological 
models/codes, e.g. the Lund model. 
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FIG. 9. Multiplicity distribution in pp collisions, at 4 
different energies from At each energy the cross section 
and multiplicity are rescaled, to put the low-N part at the 
universal KNO curve (solid line). This is done to see better 
the behavior of the "second component" discussed in the text. 

Our main suggestion for further work is that although 
the instanton-induced mechanism yields relatively small 
cross sections, this mechanism is likely to dominate over 
events with multiple color rearrangements. Is there ex- 
perimental evidence for this assertion in high-energy 
hadronic collisions? An answer is provided by Fig. ^ 
(taken from p3|] ), which shows a (specially normalized) 
compilation of multiplicity distributions in pp collisions 
at various energies. The data shows that there is indeed 
(at least) two components: (i) one, with the cross sec- 
tion <7i(s) and standard KNO distribution (well known 
from lower energy pp collisions), as indicated by the solid 
curve; and (ii) another one with a different cross section 
(72 (s) and much higher multiplicity. Ascribing the main 
peak at N / < N± >s=s 0.8 to a single color rearrangement 
reaction (=2 QCD strings formed), one can conclude that 
at the highest energy y/s — 1800GeV the multiplicity 
seen may amount up to 10 strings. 

The existence of the second component with double 
and triple multiplicity was anticipated by people doing 
Regge theory decades ago, in the form of multi-pomeron 
exchanges. However, the multiplicity data shown in Fig.|^ 
do not really fit well into this description. The second 
components simply does not look as iterations of the 
first one. There are no separate peaks and, more im- 
portantly, the s-dependence is completely different. The 
first component is in fact consistent with the approxima- 
tion used above, namely asymptotically constant cross 
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section (zero pomeron intercept), while the latter grows 
with yfs very strongly. 

Attempts to solve this puzzle in pQCD, by summing 
ladder-type diagrams in leading log(x) approximation are 
well known H, and they do indeed produce strongly 
growing cross sections and multi-parton states. So, the 
second component may well be due to those perturbative 
processes. Non-perturbative approaches (aiming mostly 
at the "soft pomeron" or the first component discussed) 
have also been tried, from old fashion multi-peripheral 
hadronic model s (e-g. recent work |34j|) to mixed gluon- 
hadron ladders |pq,pl[ . Unfortunately, none of these ap- 
proaches have lead so far to a quantitative theory. 

Results/estimates made in this work lead to the con- 
clusion, that instanton-induced color exchanges should 
dominate over pQCD t-channel gluons starting from the 
double exchange amplitudes. It is therefore logical to con- 
jecture, that the second high-multiplicity component of 
pp collisions may be generated by this mechanism. That 
would explain why multiple-string events are not just it- 
eration of the first component, and even consistent with 
where the transition appears to be. Needless to say that 
much more work is still needed for a further test of this 
conjecture. 
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